Astronomy & Astrophysics manuscript no. preprint 


February 2, 2008 


(DOI: will be inserted by hand later) 





Tracing the Sgr Stream with 2MASS 
Detection of Stream stars around Outer Halo Globular Clusters * 

M. Bellazzini^ R. Ibata^, F.R. Ferraro\ V. Testa"^ 



cn 
o 
o 

< 
oo 



> 
O 

in 

o 
m 
o 

Oh 

I ' 

o ■ 



' INAF - Osservatorio Astronomico di Bologna, Via Ranzani 1, 40127, Bologna, Italy e-mail: bellazziniObo . astro . it 

^ Observatoire de Strasbourg, 67000 Strasbourg, France e-mail: ibata@astro . u-strasbg . fr 

^ Dipartimento di Astronomia, Universita di Bologna, Via Ranzani 1, 40127, Bologna, Italy e-mail: ferraroObo . astro . it 

* INAF - Osservatorio Astronomico di Roma, Via Frascati 33, 00040, Monte Porzio Catone, Italy e-mail: 



testa@mporzio . astro . it 



Received February 25, 2003; Accepted April 24, 2003 



Abstract. We use infrared Color Magnitude Diagrams from the 2-Micron All-Sky Survey (2MASS) to search for stars be- 
longing to the tidal stream of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) around selected Galactic globular clusters. 
Statistically significant detections are presented for the cases of Pal 12 and NGC 4147, strongly supporting the idea that these 
clusters are associated with the Sgr Stream and that they were previous members of the Sgr dSph galaxy. 
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1. Introduction 

There is clear evid ence tha t the Sagittarius dwarf Spheroidal 
galaxy (Sgr dSph; llbata et al., .1994. 1997,) has lost (and is 
currently loosing) a significant fraction of its stars under 
the strain of the Milky Way tidal field. The lost stars seem 
able to remain coherently aligned with the orbital path of 
the Sg r dSph as a long - lived (> few Gyr; see Ijohnston et alJ 
( Il999l) : llbata & Lewii ( Il998l) ') stellar relic, known as the 
Sgr Stream. Parts of the Sgr Stream have been detected 
by many different groups and with many different tech- 
niques, up t o > 10() ° aw ay from the main body of the 



il fl996: Mateo. Olszg 

199i iMaiewski et all 11994 iHelmi et all Il999t llbata et al 



galaxy (Fahlman et all 11996: Mateo. Olszewski & Morrison 



2002'; Martinez-Delgado et al 
et al... 2000; Ivezjc et al 



2002'; 'Kundu et al 
2001 1 : iVivas et al., 2001; Yann 
2000l) . Direct identifications, based on the similarity of the 
CMD in the main body of the galaxy and in portions of the 
St ream ~ 160° away in g alactic longitude have been provided 

bv lNewberg et all ( 120021). 

llbata & LewisI ( Il998l hereafter IL98) and 'ibat a et all 
ll2001ah have simulated the evolution of the Sgr dSph over sev- 
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eral orbital periods (P ~ 1 Gyr), computing the orbit of the 
galaxy as well as the phase-space distribution of the debris 
under different assumptions about the flattening of the CDM 
halo. The initial conditions of the simulations were based on 
the known position and radial velocity of Sgr dSph and on its 
proper motion as estimated bv llbata et all ( Il997l l2001al) . The 
orbit has a planar rosette structure, with the pole of the orbit 
located at [f - 90°, b - -13°] (i.e. a nearly polar orbit), and 
peri- and apo-Galactic distances of 15kpc and 60kpc respec- 
tively. The derived orbit has successfully been used to predict 
the observed posi tion (Ibata et al. 2001a b) and radial velocity 
dYannv et alJ2003l) of the Sgr Stream. 

In a previous paper (Bellazzini, Ferraro & Ibata 2003 , here- 
after BFI03) we compared the IL98 orbit during the last orbital 
period with a subset of galactic globular clusters, namely those 
having lOkpc < Rqc ^ 40kpc which we termed Outer Halo 
(OH) globulars, in the phase-space of Galactocentric cartesian 
coordinates and radial velocity {X, Y, Z, V/- see BFI03 for fur- 
ther details). We found that the region of phase-space occu- 
pied by the Sgr Stream hosts a significant overdensity of OH 
clusters, thus some of them have to be physically associated 
with the Stream and were former members of the Sgr dSph. 
While the association of at least 4 OH clusters with the Stream 
is quite firmly established, the result is statistical in nature and 
we cannot confirm the actual membership of individual clusters 
with the adopted technique. Hence, we simply provided a list 
of candidate Stream members ranked according to their prox- 
imity to the Sgr orbit in the considered phase-space. Strong 
support for individual associations may be provided by a direct 
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detection of Sgr Stream stars around a candidate cluster and lo- 
cated at the same distance from us, as provided, for instance, by 
iMartinez-Delgado et al .' (2002) i n the case of our first-ranked 
candidate Pal 12 (see also .Dinescu et alJ2000l) . 

This kind of search requires the sampling of wide areas of 
sky, since the Stream is a very low surface brightness structure, 
and the use of suitable tracer stars that can be singled out from 



stars (see 


Ibata et al 


2001a; Yannvetal. 2000; Ivezic et al. 


l2000HNewbers et aL,.2002: .Vivas et al...200k Jbata et al...2002 



for the different tracers used to date). Here we use stars in 
the upper Red Giant Branch (RGB) from the near infrared 
Color-Magnitude Diagrams (CMD) of wide areas of sky pro- 
vided by the 2-Micron All-Sky Survey (2MASS) to search 
for Stream stars around the top-ranked candidates listed in 
Tab. 1 of BFI03. Sgr Stream stars at the "right" distance 
are detected around Pal 12, thus confirming the result by 
iMartinez-Delgado et al] ( l2002h . and around NGC 4147, a com- 
pletely new detection. The present study has to be considered 
as a preliminary "pilot" analysis. We will perform a deeper and 
more complete search once the final all-sky 2MASS catalogue 
is released. 

The plan of the paper is the following: in §2 we present 
the data and the method; in §3 and in §4 we describe the re- 
sults obtained for Pal 12 and NGC 4147, respectively; in §5 we 
deal briefly with the case of the other considered candidates; §6 
gives a short summary and a discussion of the reported results. 




(J-K)o 



Fig. 1. Ks,J - Ks CMD of a ~ 4° x 4° area around the center 
of the Sgr dSph galaxy. The upper RGB of Sgr stands out very 
clearly ai J - Ks > 1, to the red of all the features formed by 
galactic field stars. The lower box encloses the Sgr RGB stars, 
the upper box encloses Sgr bright AGB stars. 



2. Data and method 

All the photometric data are from the Point Source Catalogue 
(PSC) of the Second Incremental Data Release of 2MASS' 
and consist of caUbrated J,H,Ks ^ magnitudes and positions of 
sources classified as stellar and without any detected contam- 
ination from other nearby sources. Magnitudes and colors of 
all the retrieved sources have been corrected for interstellar ex- 
tin ction b y interpolating in the IRAS-DIRBE reddening maps 
of lSchlegel. Finkbeine r & Davis (1998ii) and adopting the ex- 
tinction laws bv lsavag e & Mathisl Tl979l) . The present sky cov- 
erage of 2MASS is quite patchy, hence there are regions of sky 
that cannot be scrutinized. We typically select areas of 20° x20° 
toward interesting directions, but in all cases large unobserved 
zones are included in the selected areas. Typically, the effective 
sampled area ranges from ~ i to ~ | of the selected area (see, 
for example. Fig. 4 and Fig. 5, below). Distances and other pa- 
rameters for globular clusters are taken from the compilation 
bv .Harris (1996). 

Fig. 1 displays the K,J-K CMD of a 4° x 4° field centered 
on the point of maximum density of the main body of the Sgr 
dSph. A similar field has been used by£ole QQQl) to study the 
metalUcity of the galaxy an d virtually the same field has been 
considered by'Xlara (2001) in his analysis of the metallicity 
gradient. As noted by these authors the brighter part of the RGB 

' see Cutri, R.M., et al., 2000, Explanatory Supplement 
to the 2MASS Second Incremental Data Release, 
[http://www.ipac.caltech.edu/2mass/releases/second/doc/expsup.html 

^ From now on we drop the S index from , for brevity. 



of the Sgr dSph stands out very clearly in this diagram, to the 
red of any feature due to ga lactic field s tars i n the region ( J - 
K)o > 0.8 and Kq < 12 (see^ ColellEoOll) and lAtodlll200ll) for 
discussion of the field features in the CMD). The lower part 
of the reported box encloses the more clearly visible part of 
the Sgr RGB. The upper (smaller) sector of the box encloses 
bright Asymptotic Giant Branch (AGB) stars of Sgr We will 
use only RGB stars as tracers of the Sgr stellar population but 
we maintain also the upper sector of the box to check also if 
rarer AGB stars fall in our selection box. The re ddest Sgr stars 
from 2MASS have been already used bv llbata et a l. (2002) for 
a pole-count analysis that detected an overdensity of such stars 
in the plane of the sky in the directions predicted by the IL98 
orbit. 

The detection method is quite simple: since Sgr RGB stars 
stand out clearly against galactic field stars in the consid- 
ered CMD, we will look if similar stars can be found also in 
wide fields surrounding OH globulars that are candidate stream 
members according to BFI03. The RGB feature, if present, has 
to be shifted in magnitude by the difference in true distance 
moduli between the Sgr dSph and the considered cluster (A/io). 
Therefore, we will simply report the Selection Box (SB) shown 
in Fig. 1 in the CMDs of the selected area and of a suitable 
Control Field (CF) applying a shift by A/iQ. Then we will check 
if there is a significant overdensity of stars in the SB of the se- 
lected area with respect to the CF. It is important to note that the 
SB has been constructed to include also the bluest stars in the 
Sgr RGB, in order to remain sensitive to the searched feature 
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notwithstanding the effec t of the meta llicity/population radial 
gradient observed in Sgr (lAlardil200ll) . 

Since we are looking for very low-surface-brightness struc- 
tures, it is of the utmost importance to avoid any possible 
source of contamination. To this end, within the SB we will 
select only stars having (7 - K)o > 0.95, to avoid the in- 
clusion of the reddest stars from the local M dwarf plume at 
(J - K)o ~ 0.8, and we will reject stars with Kq > 12, since 
in this range of magnitude a sparse cloud of red sources (possi- 
bly misclassified galaxies-*) provides additional contamination. 
For this reason, the adopted technique is nearly insensitive to 
the feature of interest for AfiQ i -i- 0.5, i.e. the Sgr Stream RGB 
can be detected only up to a distance of ~ 32 Kpc from us. In 
some cases (e.g., NGC 4147) the considered clusters fall in an 
unobserved zone of the selected area, thus contamination of the 
SB by cluster stars cannot occur. In all other considered cases 
we checked the CMD of the stars in the immediate vicinity of 
the cluster (up to 1° from its center) and we never found the 
case of a possible cluster star falling in the SB. This is due 
to the fact that the considered clusters are signifi cantly mor e 
metal poor th an the main population of Sgr (see Cole 2001; 
Monaco et alJl200 2 and references therein), hence their RGBs 
are much bluer and fall within the Galactic field features of the 
CMD. 

In this first attempt we considered the first six clusters of 
the ranked list of candidates of BFI03, these are, in order of 
increasing phase-space distance from the IL98 orbit: Pal 12, 
NGC 4147, Pal 5, NGC 5634, NGC 5053 and Ter 3. Searching 
the PSC catalogue we have found suitable selected areas for 
Pal 12 and NGC 4147. The present sky coverage of 2MASS 
around NGC 5634, Pal 5 and NGC 5053 is insuflicient to ob- 
tain an effective sampled area useful for the present purpose. 
Finally, we recently realized (Ma rsakov, private communica- 
tion) that the 2003 version of the iHarrisI ( 11996 ) catalogue re- 
ports significantly revised estimates of the reddening and dis- 
tance modulus of Ter 3 [E(B-V)=0.72; (m - M)v = 16.61] 
with respect to the 1996 version that was used by BFI03 [E(B- 
V)=0.32; (m-M)v = 18.00]. The new distance estimate places 
the cluster quite far from the Sgr orbit. As a consequence, Ter 3 
was removed from the list of candidates and it is not consid- 
ered in the present analysis. Thus, the test remains to be done 
for Pal 5, NGC 5634 and NGC 5053 (these cases are shortly 
discussed in §5, below), while in the following two sections we 
show the results we have obtained for the cases of Pal 12 and 
NGC 4147. 

3. Pal 12 

To search for Sgr Stream stars around this cluster, whose galac- 
tic coordinates are (I; b) - (30.512°; -47.681°), we selected the 
area comprised in the ranges 20° < / < 40° and -57° < b < 
-37°. The effectively sampled field is 71% of the selected area. 
In the upper panels of Fig. 2 the CMDs of the Pal 12 field is 
compared to the CMD of a suitable nearby Control Field at sim- 
ilar galactic latitude (80° < / < 100° and -55° <b < -35°). 

' The feature is more evident in less crowded fields with respect to 
that shown in Fig. 1; see Fig. 2 and Fig. 3. 



The difiference in true distance modulus between Pal 12 and 
the core of the Sgr galaxy is A//o = iiq{P\2) — jUo(M54) - 
-0.76 + 0.15. Hence the SB shown in Fig. 2 has been shifted 
by AyUo - -0.76 in K to search for the RGB feature at the dis- 
tance of Pal 12. 

Even from a visual inspection it is clear that Sgr RGB stars 
are present in the field suiTounding the cluster and are miss- 
ing in the CF. To quantify the observed difference we count 
the stars falling into the SB (with the additional conditions 
(J - K)o > 0.95 and Kq < 12.0) in the two CMDs (Nf" 
and A^^^, for the cluster field and for the CF, respectively). 
To take simultaneously into account the differences in sam- 
pled area and any possible difference due to galactic popula- 
tion gradients between the cluster fields and the CFs we nor- 
malize to the number of stars spanning the same magnitude 
range of the plotted SB and with colors 0.55 <( J - K)o < 0.90 
(^iwrm ^iid ^norm)- Notc, howcvcr, that all the results are fully 
confirmed if we normalize by the ratio of the sampled areas 
instead. In particular, the normalization ratios by number and 
by area are virtually identical in all of the considered cases. 
Defined R - N'^nmml^norm^ computc the normaUzed differ- 
ence: 

AA^* = Nl'"' - X Af^ , (1) 

and its standard deviation, assuming Poisson eiTors in star 
counts: 

cr(AA*) = ^Nl'"' +R^x A^^ . (2) 

In the present case AA* - 53.5 + 10.9 is obtained, that is, a 
very significant diff'erence (4.9cr). The searched RGB feature is 
indeed present around Pal 12, it is statistically significant and it 
is found at the right magnitude level. 

To explore better the nature of the detected signal we repeat 
the counts shifting the SB in the range AyUo = -0.76 + 1 .0 at 
0.04 mag steps. We plot the resulting AA* as a function of Ajuo 
in the lower panel of Fig. 2. The normalized difference reaches 
a wide maximum in the range -1.1 < A/zq ^ -0.5 and the 
characteristic Afio of Pal 12 (indicated by a filled triangle in the 
top of the panel) falls right in the middle of the detected plateau. 
The plot has to be considered just a consistency test since the 
constraining power on distance of a sparsely populated upper 
RGB sequence is obviously not strong. Moreover, the signal 
remains significant over a wide range of AA* because, after 
all, the Sgr stars are present in the Pal 12 field and absent in 
the CF in this area of the CMD. Thus any position of the SB 
that includes some of the Sgr RGB stars provides a sensible 
difference with respect to the CF. The effects of this syndrome 
may be particularly strong in shifting the SB toward brighter 
magnitudes (Kq ~ 10.5) since in this region the CF is virtually 
devoid of stars. 

Notwithstanding the above caveats the detection of Sgr 
Stream stars around Pal 12 is clear and signifi cant, confirming 
the findings by Martinez-Delgado et al ] ll2002h . Given it's posi- 
tion in the phase-space with respect to the orbit of the Sgr dSph 
(see BFI03) and its physical location within the Sgr Stream, we 
conclude that Pal 12 presently belongs to the Sgr Stream and 
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Fig. 2. Upper panels: the CMD of the selected area around the position of Pal 12 (left upper panel) and of a Control Field located 
at similar Galactic latitude. The reported Selection Box has been shifted by A/uq = -0.76 with respect to Fig. 1, to displace it to 
the distance of Pal 12. The selected stars are marked with heavier points. Lower panel: The normalized difference in the number 
of stars falling in the Selection Box as a function of A/io- Different symbols correspond to different values of AA^* in units of cr 
(filled points: AN-k/cr < 3 , triangles: 3 < A^*/cr < 4, squares: 4 < A^*/cr < 5, etc.). The large filled triangle on the top of the 
panel marks the assumed A/io for Pal 12. 



that it was a former member of the globular cluster system of 
the Sgr dSph. 

Even an educated guess of the surface brightness (S) of 
the Sgr Stream at the considered locations would be of great 
interest for future modeling of the dynamical history of the 
dwarf galaxy. We obtain such an estimate by counting the num- 
ber of stars in the SB having luminosities within 1 magni- 



tude of the assumed RGB Tip (N^j'^, dividing by the sam- 
pled area and rescaling by the central surface brightness of Sgr 
dSph to the ratio of surface density of RGB Tip stars in the 
considered area and in the center of Sgr The underlying as- 
sumption is that the stellar population in the center of the Sgr 
dSph and in the Sgr Stream are similar. The adopted method is 
widely used (see, for example. . Klevna et al.. 1998: Mateo et aL 
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Fig. 3. The same as Fig. 2 for the selected area around the position of NGC 4147 and of a Control Field symmetrical with respect 
to the Galactic equator. AyUo - -0.70, in this case. 



Il99ll) and its physical basis is anchored on the Evolutionary 
Flux Theo rem {Renzini & Buzzoni 1986; Renzini & Fusi Pecci 
ll988tlRenzinill998l) . In the central 0.5 x 0.5 deg^ of Sgr dSph 
we find Nj.'^ - 50 stars within 1 mag from the RGB Tip and 
we assume that this corresponds to the central surface bright- 
ness o f £v(0) ~ 25.4 + 0.3 mag/arcsec^ reported bv iMateol 
lll998l) . The possible contamination by field stars is tentatively 
removed with the same technique used above, i.e. considering 
nI''' = A^^'^^' -Rx N^[^ instead of N'^'^'^ itself. Finally, we ob- 
tain for the average surface brightness in the Sgr Stream around 
Pal 12, Ey - 33.9 + 0.5 mag/arcsec^. The uncertainty has been 



estimated by 1000 Montecarlo simulations in which lly^'(0), 
Nf;^^ and N^f^ were simultaneously extracted from gaussian 
distributions having the observed values as mean and the re- 
ported 1 - cr uncertainty as cr. 

4. NGC 4147 

NGC 4147 [il; b) = (252.850°; 77.189°)] Ues in a region of the 
sky (220° i I i 320°) where the orbit of the Sgr dSph runs ap- 
proximately at constant galactic latitude (b ~ 75°). We take 
advantage from this occurrence by taking a wide selected area 
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comprised within 230° < I < 290° and 65° < b < 85°. 
The effectively sampled field is 50% of the selected area. As 
Control Field we choose the region of the sky symmetric with 
respect to the galactic equator, i.e. 230° < / < 290° and 
-85° < b < -65°. The resulting CMDs are shown in Fig. 3, 
which is arranged in the same way and with the same sym- 
bols as Fig. 2, above. The diff'erence in true distance mod- 
ulus between NGC 4147 and the core of the Sgr galaxy is 
Ajuo = //o(A^4147) -//o(M54) = -0.70 + 0.15 and the SB has 
been shifted accordingly in these panels with respect to Fig. 1. 

Also in this case the Sgr RGB feature is quite evident and 
significant. The normalized difference is A^V* - 35.4 + 8.7, a 
4. Oct signal. The lower panel of Fig. 3 shows a less pronounced 
maximum with respect to the Pal 12 case, probably due to the 
higher degree of contamination of the considered CMDs in the 
region (J - K)o > 0.95. The adopted AjUq for NGC 4147 falls 
at the "bright" edge of the plateau. However if we consider the 
normalized difference in terms of standard deviations {N^jcr), 
the maximum signal is detected at AyUo = -0.70 + 0.04, i.e. at 
the distance of NGC 4147. 

Hence, it has to be concluded that also this cluster is im- 
mersed in the Sgr Stream and that it is physically associated 
with this structure. 

Using the technique illustrated above we obtained = 
34.0 + 0.6 mag/arcsec^ for the Sgr Stream surrounding 
NGC 4147. The comparison with the analogous estimate ob- 
tained in the region around Pal 12 suggests that the Stream 
maintains nearly c onstant surface brightness along its whole 
extension (see also lMateo. Olszewski & Morrisonll998i) . 

5. Other cases 

In this section we shortly comment on the other candidates by 
BFI03 considered by the present study. In all cases we cannot 
reach a firm conclusion, either because of a too small effec- 
tively sampled area or because of problems with interstellar 
extinction. 

NGC 5634 and Pal 5 are relatively nearby in the sky and 
have similar A/io, and thus may be studied simultaneously. 
Fig. 4 shows that the coverage of 2MASS in this region of 
sky is quite poor and unfavorably placed with respect to the 
considered clusters and to the Sgr orbit. Hence, the absence 
of a significant detection in the selected area shown in Fig. 4 
cannot be considered as evidence that the clusters are not as- 
sociated with the Sgr Stream. Hopefully, the question should 
be resolved when future releases of the 2MASS catalogue are 
available. 

The same problem plagues the application of the test to 
NGC 5053. In this case only a barely significant detection is 
obtained but the effective sampling is clearly insufficient and 
unfavorable (see Fig. 5). 

In Fig. 6 we report, for comparison, the distribution of the 
sampled fields in the selected area around NGC 4147. 

6. Summary and Conclusions 

We have examined the 2MASS Kiy,{J - K)q CMDs of wide 
areas of the sky around six Outer Halo globular clusters to 




320 340 360 



Fig. 4. The fields sampled by 2MASS in the selected area 
around NGC 5634 and Pal 5 (in grey). The continuous line is 
the projection in the sky of the Sgr orbit by IL98. 
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300 320 340 360 



Fig. 5. The fields sampled by 2MASS in the selected area 
around NGC 5053 (in grey). The continuous line is the pro- 
jection in the sky of the Sgr orbit by IL98. 

search for stars belonging to the Sgr Stream and that are lo- 
cated approximately at the same distance as the considered 
clusters. The considered clusters were selected for their phase- 
space proximity to the orbit of the Sgr galaxy, according to 
BFI03. The dataset was found to be suitable for the apphca- 
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Fig. 6. The fields sampled by 2MASS in the selected area 
around NGC 4147 (in grey). The continuous line is the pro- 
jection in the sky of the Sgr orbit by IL98. 

tion of the test to the cases of Pal 12 and NGC 4147. In both 
cases, an overdensity of tracer stars of similar heliocentric dis- 
tance as the clusters, was detected at the > 4.5cr level of sig- 
nificance with respect to Control Fields. Hence we have shown 
that both clusters are physically immersed in the Sgr Stream 
and are very likely members of this structure as well as be- 
ing former members of the Sgr dSph galaxy. While our re- 
sult for Pal 12 confirms a previous (independent) detection by 
[Martinez-Delaado et al. (2002), this is the first time that Sgr 
Stream stars are found around NGC 4147. A thorough analysis 
of the structure and the stellar populations of the Sgr Stream 
will be presented when the full-sky 2MASS catalogs become 
available. 
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